Evading l/m;,-suppressed IR divergencies in QCDF: 
Bs KK Decays and Bd^s mixing* 
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We analyze the deviations of the mixing induced CP asymmetry in _B (j>Ks from sin 2/3, as well 
as the deviations of the asymmetries in Bs K*^K*'^ , Bs 4>K*° and Bs 4>(j) from sin2/3s, 
that arise in SM due to penguin pollution. We use a theoretical input which is short-distance 
dominated in QCD-factorization and thus free of IR-divergencies. We also provide alternative ways 
to extract angles of the unitarity triangle from penguin-mediated decays, and give predictions for 
Bs -» K*°K'-° observables. 



I. INTRODUCTION 

The phenomenology of hadronic i?d-decays has been 
a matter of intensive research in the past 15 years, in 
part due to the large amount of data collected at the 
B-factories Babar and Belle, and at CDF and DO. This 
research has led to the consensus that the CKM mech- 
anism for CP and flavor violation is accurate. However, 
some puzzles have survived up to this day [l|, H, H, 0, IB ■ 

A new era in B-physics has been triggered by the ex- 
perimental precision and by the emergent exploration of 
the Bg system. The future is marked by the starting 
of LHCb and the possibility of a super-B factory. So 
far, measurements on the Bs system are limited to the 
mass difference AMg [6l|, other mixing parameters ^j, 
and some Bs — > ttK and Bg — > KK modes but the 
extension of this list is an important element in the B- 
physics program. 

On the theoretical side, the study of non-leptonic B- 
decays is difficult because of the presence of important 
long distance strong interaction effects. The correct com- 
putation of these contributions is crucial to be able to 
resolve small NP contributions. The amplitude of a B 
meson decaying into two light mesons can be written as 



A{B -> AhM2) = Ai^)* Tm,m, + A(^)* 



Ml Ma 



(1) 



were Ag^^ = VqbV*jj, and T and P are called "tree" and 
"penguin" . These hadronic parameters can be extracted 
from data, or can be predicted using symmetries (such 
as flavor). The direct computations from QCD are much 
more involved and are based on factorization and the 
l/rrih expansion. They appear in the context of QCDF, 
pQCD or SCET @. While the methods based on fla- 
vor symmetries include naturally all kinds of long dis- 
tance physics, they suffer from big uncertainties due to 
bad data an poorly estimated SU(3) breaking. On the 
other hand, methods based on factorization suffer from 
uncertainties due to non-factorizable chirally enhanced 
l/rrib corrections, and long distance charm loops (charm- 
ing penguins). 
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These proceedings review a recent proposal to improve 
(at a phenomenological level) on some of the weak points 
of the approaches mentioned above [l^, [U- We also 
include a straightforward application to -B ^ (t>Ks and 
comment on the limitations and the applicability of the 
method. 



II. EXPRESS REVIEW OF Bq ~ Bq MIXING 

The time evolution of a Bq meson {q — d, s) can be 
easily described by changing to the mass eigenbasis. The 
relationship between the flavor basis and the physical 
mesons is specified by a mixing parameter usually de- 
noted by q/p, 



\Br) = 



l + la/pPV 



VI + \1/P\ 
\Bh) = ^=1==(\B-)~:^\B 
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(2) 



The time evolution of the mass eigenstates is straight- 
forward and depends only on the masses and the widths 
of the physical B-mesons. Therefore, the evolution of 
the flavor eigenstates -which describes the oscillations- is 
specified by the masses and widths of the physical mesons 
and on the mixing parameter q/p. 

Then one can define the mixing angle 4>m as 

0M = - arg (q/p) (3) 

In terms of the entries of the effective hamiltonian, 
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where we have used that |ri2| ^ |Mi2|. This means that 
the above definition of the mixing angle is equivalent to 



arg (M12) 



(5) 



One should be aware, though, that these quantities are 
not convention independent and are sensitive to unphys- 
ical phase redefinitions. However, once a convention 
for the weak phases is chosen everywhere, this defini- 
tion is meaningful. In the Wolfenstein parametrization 
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Vub and Vtd have phases of 0(1) and Vts of O(A^). 
This means that in SM, in the case of Bd ~ Bd mix- 
ing oc (V;*Ftfc)2 and (f>j^'' = 2/3 + 0{X^), and in the 
case of Bs - B^ mixing Mf^ oc (V^* t/tb)^ and (f>^^'' = 2/3s. 

In order to measure the mixing angle one usuaUy 
looks at CP asymmetries. The time-dependent CP asym- 
metry oi a B{t) ^ f decay is defined as 



T{B{t) ^ /) - T{B{t) ^ /) 
riB{t)^ f)+TiB{t)^ f) 



(6) 



where B{t) {B{t)) was a B (B) aX t = 0. For / = fcp 
a final CP-eigenstate (with CP eigenvalue 77/), and ne- 
glecting the small CP violation in mixing one finds 



Acpit) = 



Adh- cos(AAft) + Ami^ sm{AMt) 
cosh(Art/2) - Aaf sinh(Art/2) ' 



(7) 



where the quantities Adh and -4inix are the so called di- 
rect and mixing induced CP asymmetries, and are mea- 
sured from the time oscillations of Acp{t)- The mixing 
induced CP asymmetry is given by 
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where A/ is a convention-independent (physical) quan- 
tity. Here Aj = A{B ^ f) and Aj = A{B -> /). 

In some particular cases, one can extract the mix- 
ing angle in a very clean way from a mixing induced 
CP asymmetry. The prominent example is the case of 
Bd J/tpKs Since this decay is dominated by a 

single amplitude, Aj/^kJAj/^k^ ~ Vj/i^K, = -1, and 
therefore, 



- AmUBd -> J/^Ks) ~ sin 2^ 



(9) 



The neglected amplitude is both CKM and as{mb) sup- 
pressed with respect to the dominant amplitude, so the 
corrections to this equation are below the percent (or 
even the per mil) level [Ulillis'l. 

The same argument is valid in principle for penguin- 
mediated b —^ s decays, where the amplitude can be writ- 
ten as 



A = A(f)* T + Xi'''* P ~ a1")* p. 



(10) 



The tree part is strongly CKM suppressed since 
|Au^/Ai^^| ~ 2%. However, as opposed to the previous 
case, T is not suppressed with respect to P. In fact, an 
enhancement of T over P would spoil the extraction of 
4>M from these modes. 

The classical example is the case of Bd <j>Ks , where 



In order to have under control the smallness of the cor- 
rection AS^Ks o'^s should be able to bound the size of 
^e{Tij,KjP4,Kj- In fact, latest data gives jl6| 



AS-^^^ = -0.39 ± 0.20 



(13) 



so if the uncertainties are reduced around this central 
value, the claim of a NP signal will have to rely on a 
solid bound for the SM penguin pollution. 

A first approach to bound this tree-to-penguin ratio 
was taken in the context of flavour SU(3) symmetry 
[TtI [Tsj . The argument is that if there was a large hi- 
erarchy between T and P in a 6 s mode, the hierarchy 
would persist when one moves to the SU(3)-related 6 — > d 
modes. But in these modes the tree is not suppressed 

with respect to the penguin, because |Ai'*''| |Ai'^''|, so 
this hierarchy would have a clear impact in the observ- 
ables of the b ^ d modes. In this way one can write the 
following bound 



\AS, 



4>K^ 



< V2X 




valid in the SU(3) limit and under a non-cancelation as- 
sumption between Bd <j>Ks and (t>K~^. With 
the present data the bound is 



|A5|^f)|<0.4 



(15) 



The same kind of analysis can be applied to other 
penguin-dominated decays [H, [l^ . 

A second approach has been followed in the framework 
of QCD-factorization ^2d\ , which gives a much more com- 
petitive bound. 



0.01 < AS^™^ < 0.05 



(16) 



Related analyses have been carried out in SCET [2l| and 
pQCD [131 . For a recent review see also [23j . 

Concerning the Bs — Bs mixing angle, the clean tree- 
level determination comes from Bg J/ipcf). The re- 
lated penguin-mediated decays are, for example, Bg — > 
K*°K*°, Bs <j>K*° or Bs (jx/). Again, one can write 



VfAmUBs f) = sin 2/3, + ASj. 



(17) 



A study of the amounts by which their mixing induced 
CP asymmetries deviate from sin2/3s in the SM can be 
found in [111,!!. 

In the next pages we follow the approach in [l^, [IH , 
based on a theoretical input that we call A. 



III. THEORETICAL INPUT 



{Bd ^(pKs) = sin 2/3 + AS", 



with 



A^, 



As) 



Re 



p 



<I>K, (11) Consider the quantity A = T — P. This quan- 

tity is a hadronic, process-dependent, intrinsically non- 
perturbative object, and thus difficult to compute theo- 
retically. Such hadronic quantities are usually either ex- 
sin 7 cos 2(3 + ■ ■ ■ (12) tracted from data or computed using some factorization- 
based approach. In the latter case, A could suffer from 
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(2.29 ± 0.67) X 10 


GeV 




(1.85 ±0.79) X 10" 


-"GeV 




(1.62 ± 0.69) X 10" 


""GeV 




(1.16 ± 1.05) X 10" 


-"GeV 




(2.06 ± 2.24) X 10" 


-"GeV 



TABLE I: Values of A for the various decays of interest. In 
the case of two vector mesons these numbers correspond to 
longitudinal polarizations. 



the usual problems related to the factorization ansatz 
and in particular long-distance effects. 

However, for a certain class of decays, T and P share 
the same long-distance dynamics: the difference comes 
from the (u or c) quark running in the loop, which is 
dominated by short-distance physics ■ Indeed, in such 
decays, A = T — P is not affected by the breakdown of 
factorization that affects annihilation and hard-spectator 
contributions, and it can be computed in a well-controlled 
way leading to safer predictions and smaller uncertain- 
ties. 

Table [J shows the values of A for our cases of inter- 
est. This quantity was used to predict branching ratios 
and asymmetries in KK modes, and the out- 

come was promising [13, HI]. In |26| this quantity was 
used to extract the angle a of the unitarity triangle from 
Bd K^K^ . In the following section we review the for- 
mulae that allow to extract T and P from data and the 
theoretical input A. 



IV. TREE AND PENGUIN CONTRIBUTIONS 

We begin writing the two self-conjugated amplitudes 
in terms of tree and penguin contributions, 

A^XW*rj.^^(D)*p ^ A^X(f)T + Xf^P (18) 
Now we put T = P — A and we square the amplitudes. 
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FIG. 1: Allowed values for the longitudinal direct CP asym- 
metry of Bd — > K'^K*^ in terms of its longitudinal branching 
ratio, according to the value of /S.k'k* ■ 



for two non-identical particles in the final state. The 
resulting expressions are 



PP(l+^dir)/5PS 
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BR{1 - AA,,)/gps 
|Af)+Al^^P 
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(22) 



These are the equations for two circles in the complex P 
plane, whose solutions are the two points of intersection. 
This will result in a two-fold ambiguity in the determi- 
nation of P (and T). Before writing down the analytical 
solutions, notice that in order for solutions to exist, the 
separation between the centers of these circles must be 
smaller than the sum of the radii but bigger than the 
difference. This translates into a consistency condition 
between BR, Ad„ and A: 



l^dirl < 



2BR 



2- 



7^f,A2 



2BR 



(23) 



= |A(^'* + Ai^'*|2 
\A\- = |A(^)+Al^)p 



P 



P 



a<°'+a!." 



But the squared amplitudes are directly related to ob- 
servables. 



where BR = BR/gps and TZjj is a specific combination 
of CKM elements (see Table HIl . This condition turns 
(19) out to be highly nontrivial. For example, FigUshows the 
allowed values for the longitudinal direct CP asymmetry 
of Bd 7^*0^*0 jj-^ terms of its longitudinal branching 
ratio. It can be seen that for BR ^ 3 x 10~^ the direct 
CP asymmetry is quite constrained. 

The hadronic quantities P and T arc then given by 



= BR{l+Adn.)/gps 
= BR{l-Adn)/gps 



(20) -^"^[^] = 



where gps is the usual phase-space factor. Neglecting the 
masses of the light mesons with respect to the B mesons, 



gps{Bd) = 8.8x lO^GcV-^ 
gps{Bs) = 8.2 X lO^GcV"^ 



(21) 



BR Adir 



T = P + A 



\ 



-/m[P]2 - 



BR 



(24) 



4 



Co 


„(d) 
^1 


^2 


7^d 


-3.15 • lO"'^ 


-0.034 


6.93 ■ IQ-^ 


7.58- 10"^ 


J") 
Cq 


''I 


^2 




3.11 ■ 10"^ 


0.011 


1.63 • 10"^ 


1.54 • 10"^ 



TABLE 11: Numerical values for the coefficients c^^' and TZd 
for 7 = 62°. 



The decay Bg x*'^K*'^ is more appropriate because 
^K-'K' is under a much better numerical control, but 
there is no experimental value for the branching ratio. 
For the sake of illustration, we just mention that for 



(30- 40) X 10^^ one gets 



0.037 < ASk'K' < 0.051 



(30) 



where the coefficients c|^^ are again some specific com- 
binations of CKM elements (see Table HI]) . 

Equations ([24|) allow to extract the hadronic parame- 
ters T and P from experimental data on BR and AdiT, in- 
formation on sides of the unitarity triangle and the weak 
phase 7, and the theoretical value for A. This method is 
also powerful because if no experimental information is 
available for ^dir, one can just vary ^dir over its allowed 
range in eg. ((23)) . So in fact T and P can be extracted 
from BR, A and CKM elements. 



VII. OTHER ANGLES FROM DATA AND A 

There are also some expressions that can be written 
that relate directly branching ratios and asymmetries to 
other angles of the unitarity triangle through the quan- 
tity A JJjI. These expressions do not require any CKM 
angle as an input, just sides of the unitarity triangle. 
The experimental input is minimized by measuring , 
which can be extracted from the time-dependent un- 
tagged rate [H, |2^. Then, in the case of a Bd meson 
decaying through a b ^ D process {D = d,s), 



V. sin 2/3 FROM B -» (j>K, 

Following the discussion in the previous section, the 
bounds for Tie{T^Kj P4>Ks) are given by 



Re 
Re 



(^) < l+{^c['^+C{BR,K.,Ai^jy' 



> 1 



C{BR^K^,AiKj 



C{BR,A) ^ v^-(4^V4'^)^ + (l/4^))5i?/A2 (25) 



Introducing the numbers for the coefficients from Table 

4>K, 



im the value or A^^ in Tableland the latest experi- 



mental value for the branching ratio [16] 

BR{Bd ^ 0if.)exp = ^■it\l X 10-6 
we get the following bounds for AS^k^ , 
0.03 < AS^K, < 0.06 



(26) 
(27) 



VI. sin 2f3s FROM Bs 



VV 



Equations (|25p apply as well to Bs ~> VV. Here we fo- 
cus on longitudinal polarizations for which the numerical 
values of A are under control. As mentioned above, pen- 
guin mediated Bg — > VV decays measure sin2/3s, but no 
experimental information is available yet for CP asym- 
metries in Bs decays. There is, though, an experimental 
number for the Bg — + (fxj) branching ratio (l6j , 



BRiBs 



)cxp = 14t? X 10-6 



(28) 



If we suppose that the longitudinal polarization fraction 
is ft"^ - 50% as QCDF suggests jl^, then we find 



0.006 < ASm < 0.072 



(29) 



sin^ a 



BR{1 - Aat) . „.„2 a _ BR{l^_A^ 



2|A 



(■D)|2 



A\ 



sin^ /3 = 



2|A^ 



(31) 



and in the case of a Bg meson decaying through a,h —t D 
process {D = d,s), 



sin Ps = 



BR{1 - Aaf) 



2IA 



(32) 



VIII. PREDICTIONS FOR Bs 



K'K* 



Assuming no new physics contributing to Bd 
K*'^K*'^, we can also give SM predictions for the branch- 
ing ratios and CP asymmetries of its U-spin partner 
Bs K*°K*" ll]. The inputs are A^.^., SU(3) 
breaking from QCDF, and the experimental value for 
BR{Bd K*°K*°). While there is still no experimental 
information on this branching ratio, it is remarkable that 
for BR{Bd K*"K*") > 5 x 10"'^ the results are almost 
insensitive to it. Taking 7 = (62 ± 6)° and 2/3^ = -2°, 
we find 

/Bi?i°"s(Bs ^ K*°K*") 



= 17 ±6 (33) 
K*"K*")sM = 0.000 ± 0.014 (34) 
K*^K*°)sM = 0.004 ± 0.018. (35) 



IX. CONCLUSIONS 

To conclude, we would like to comment on the rele- 
vance of the proposal. First, the applicability of the ap- 
proach has to be checked individually for each mode: it 
can only be applied to those decays for which A receives 
no contributions from annihilation or hard spectator 
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scattering graphs, such as S — > K^*^K^*\ Bd — > cf)K^*\ 
_B+ — !■ 7r+0, etc. The predictions derived in this way 
include most of the long distance physics, which is con- 
tained inside the experimental input. The used theo- 
retical input is minimal, and is the most reliable input 
that QCDF can offer, free from the troublesome IR- 
divergencies. Moreover, the theoretical error is under 
control and is likely to be reduced in the near future due 
to, for example, the fast progress that is taking place in 
lattice simulations. 

There is, however, an honest criticism due to the fact 
that some long distance effects that are controversial in 
QCDF, are also absent here. The prominent one is the 
contribution from the charming penguins, which could 
easily account for a significant discrepancy between the- 



ory and experiment that would not be due to NP. A 
recent example of this is given in Q. 
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